Are there sharp fractional charges in Luttinger liquids? 
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We examine charge fractionalization by chiral separation in a one-dimensional fermion system de- 
scribed by Luttinger liquid theory. The focus is on the question of whether the fractional charges are 
quantum mechanically sharp, and in the analysis we make a distinction between the global charge, 
which is restricted by boundary conditions, and the local charge where a background contribution 
is subtracted. We show, by way of examples, that fractional charges of arbitrary values, all which 
are quantum mechanically sharp, can be introduced by different initial conditions. Since the system 
is gapless, excitations of arbitrary low frequency contribute to the fluctuations, it is important to 
make a precise definition of sharp charges, and this we we do by subtraction of the ground state 
contribution. We very briefly comment on the relevance of our analysis for proposed experiments. 

PACS numbers: 



I. INTRODUCTION 

Charge quantization is regarded as a fundamental phenomenon of nature, and it has an overwhelmingly strong 
empirical support. All elementary particles that have been observed thus come with a charge that is a multiple of 
the fundamental unit charge. Although charge quantization is generally considered to be a fundamental quantum 
mechanical effect, there is no unique way to derive it from fundamental principles. There are, however, certain general 
arguments for such a connection. One set of arguments is based on the (hypothetical) existence of magnetic monopoles. 
A consistent quantum theory of systems that contain particles with electric and with magnetic charges requires the 
Dirac quantization condition to be satisfied, and this in turn implies that both electric and magnetic charges must 
be quantized [T]. Furthermore, the natural description of the electromagnetic field in the presence of monopoles is 
in terms of topologically non-trivial fibrebundlesP]. Another set of arguments is based on the assumption that the 
gauge group of electromagnetism is not the translation group M but rather the unitary group U{1), which is compact. 
This means that the electromagnetic potentials are angular variables, and from that follows that the electric charges 
are integer multiples of a fundamental unit of charge just as a quantum mechanical angular momentum is quantized 
in units of h/2 [5]. This is a natural scenario in certain unified theories where the gauge group of electromagnetism is 
embedded in a lager non-Abelian symmetry group. It is interesting that topology is a basic ingredient in both these 
approaches, and that they are closely connected in that many unified gauge theories have soliton solutions describing 
magnetic monopoles HJ |S] . 

In spite of the fundamental character of the charge quantization condition, charge fractionalization is a meaningful 
and interesting phenomenon that may take place under certain conditions, typically on the background of a non- 
trivial quantum many-particle state. It is then important to appreciate the distinction between fractionalization of 
the fundamental charge and fractionalization of the expectation value of the charge, which may arise from quantum 
fluctuations in the charge density. An example of the latter is the charge distribution of a single electron in the 
presence of two clearly separated positive ions. While the the expectation value of the local charge close to one 
of the ions will be e/2, the fluctuations are large and any local measurement of the charge will yield either or 
e. Charge fractionalization on the other hand, amount to having sharp local charges, meaning that a suitable local 
charge measurement will yield e/2. 

There are two cases of charge fractionalization in condensed matter systems that has caught much attention. One 
is the case of half-integer charges in certain one-dimensional crystals characterized by a Peierls instability. In a field 
theoretical description of this system, the fractional charges are associated with soliton excitations [SI 13 IE]- ED] The 
second case has to do with quasi-particle excitations in the two-dimensional electron gas of the (fractional) quantum 
Hall effect. In this case the excitations appear with fractional charges as well as fractional statistics, and the values 
they take are determined by the ground state, with different rational values for each plateau of the Hall conductivity 

[niiniiii]. 

An important property of the ground states where fractionally charged excitations have been predicted to exist is 
that there is a gap in the energy spectrum. For the plateau states of the quantum Hall effect this amounts to the 
two-dimensional electron gas behaving as an incompressible liquid. The effect of the gap is to suppress low-frequency 
fluctuations, and this is what makes it possible for the excitations to appear with sharply defined particle properties 
in the quantum mechanical sense |12|, 113] . 
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It is against this background, one should view recent theoretical claims of the existence of fractional charges in 
one-dimensional fermion systems described by Luttinger liquid theory, such as quantum wires or carbon nano-tubes. 
In particular, it has been pointed out that fermions that are scattered on impurities [H] or simply injected into the 
one-dimensional system [15' will dynamically be separated into a right-handed and a left-handed component, where 
each of these carry a non-integer fermion number. It has also been argued that these charges are quantum mechanically 
sharp observables |15j . which indicate that one could view them as carried by quasi-particlcs with fractional charge 
(and statistics). There have also been suggestions of how one should experimentally detect such fractionally charged 
objects ^171 US]. However, since these systems are gapless and also since the charges depend continuously on an 
interaction parameter, g, the question of whether the fractional charges can be viewed as quantum mechanical sharp 
is less clear. In particular the sharpness of the charges seems not protected by any topological argument. 

These questions are far from settled, and the motivation for the present paper is to examine more carefully in what 
sense fractional fermion numbers, created by chiral separation, can be viewed as being quantum mechanically sharp. 
We analyze the question by following Heinonen and Kohn |19j in introducing a distinction between the total local and 
global charges. The global charge is the full charge of the system, including boundary charges in a one-dimensional 
system with boundaries. For a compact space, which is here a circle since we are in one dimension, it is simply the 
integrated charge, which is the q = (momentum) component of the charge density. In the bosonic field theory 
description, the quantized values of the global charges are restricted by periodicity conditions and can be viewed 
as topological quantum numbers that are independent of the chiral separation. In the fermion picture they carry 
information about the fundamental fermions of the theory and are independent of the particle interactions. The local 
charges on the other hand are insensitive to the periodicity constraints and may take values different from those 
of the fundamental fermions. They are defined as the g — > limit of the charge densities rather than the q — 
component. We have already mentioned the principal difficulties with defining sharp local charges in gapless systems. 
Nevertheless, we shall see that it is possible to consistently subtract the contribution emanating from the ground state 
fluctuations, and we shall call a charge sharp if there are no additional contributions. 

The paper is organized in the following way. We first introduce the one-dimensional model, which in the low-energy 
approximation is equivalent to the Tomonaga-Luttinger model, and summarize the basic facts about its bosonized 
form. We next discuss chiral separation of excitations with non-integer fermion numbers by three specific examples. 
The first one is the situation discussed by Pham et. al. [15 , where a fermion is injected at one of the two Fermi 
points of the system, and we show how to reproduce the fractional charge values of the chiral components in our 
formalism. In the next example we follow the approach of Heinonen and Kohn [19 , who introduce what they call the 
local quasi-particle charge by inserting a fundamental particle and then adiabatically turning on the interaction. The 
third example describes a situation where a local polarization charge is introduced by applying an external potential. 
We show that in all these three cases the chiral components carry fractional charge, but the values differ between the 
cases. In the first two examples it depends (in a different way) on the coupling strength, while in the last example it 
also depends on the strengths of the external potential. 

Next we consider the interesting question whether the fractional charges can be considered as sharp. This we do 
by evaluating the variance of the local charge in the bosonic description, and we focus first on the ground state 
fiuctuations for which we derive an explicit result. The charge fiuctuations of the excited states are then evaluated for 
each of the three examples, and we find that in all three cases the charges are sharp in the meaning that the variance 
of the charge operator is identical to the ground state variance. This cast some doubt on the conclusion that any of 
these should be seen as the charge of a basic charge-carrying object of the Luttinger liquid. We conclude with some 
comments about this and also the possible relevance of our results for proposed experiments. In the appendix we 
compare the results of the bosonized theory with explicit calculations using free fermions and without making any 
low-energy approximations. 

Note that charge is throughout the paper taken to be dimensionless, and equal to the fermion number of the 
many-particle system. 

II. THE MODEL 

We consider a system of (non-relativistic) spinless fermions in one dimension, with a many particle Hamiltonian of 
the standard form 

H = J dx^\x)i-^^mx) + ljd^j dx'p{x)V{x - x')p{x') (1) 

where •4'{x) is the fermion field operator, p{x) = {x)ip{x) is the particle density and V{x — x') a two-particle 
interaction. The system is confined to a ring of length L, which we assume to be much longer than any relevant 
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physical length, and whenever convenient we may therefore take the limit L ^ oo. The field operator is assumed to 
be periodic on the ring, 'ip{x + L) = ip{x). 

In momentum space the Hamiltonian takes the form 



q,ki,k2 



with 

i'ix) =4fJ2 ' Ck = ^ f dx V'(a;)e-"=- (3) 



L ^ ' \/L Jo 

and where k takes the discrete values, k = 2n7r/L with n as an integer. 

With TVq particles in the system, the ground state of the non-interacting system defines a full Fermi sea, where all 
momentum states between the two Fermi points -iikp are occupied and other single particle states are empty. The 
Fermi momentum is related to the particle number hy kp = Nqtt/L. The particle interaction is assumed not to change 
this picture in an essential way. With the interaction turned on, the system in the ground state should still define a 
Fermi sea, but now with a smooth transition from full occupation to no occupation for momentum states in an interval 
Afc about the Fermi points ±kp. We assume Afc << kp and consider in the following a low energy description of the 
system which is restricted to processes that excite particles only in this interval. The particle number we also assume 
to be close to Nq and we define N to measure the particle number relative to this state. 

The low-energy approximation to the Hamiltonian takes the form of the Tomonaga-Luttinger model [501 121] 

H = vphJ2ixk- kp) : c\kCx,k ■ +^ E(^i(9)px,<7Px,-<z + ''^2(g)px,<zP-x,-<z) (4) 

x,fc x.« 

where excitations close to the two Fermi points ±kp are assigned different values for the quantum number x — il- The 
interaction is now separated in two parts, with Vi as the interaction between pairs of particles close to the same Fermi 
point, and V2 as the interaction between particles at opposite Fermi points. The form ([T|) of the original Hamiltonian 
introduces the restriction Vi = V2, however for more general (non-local) interactions the two interaction potentials 
may be different. In the low-energy approximation there is no interaction matrix element that will change the relative 
number of fermions at the two Fermi points. This means that there are two conserved fermion numbers N^, that 
both take integer values. The parameter vp in the above expression is the Fermi velocity. In the non-interacting case 
it is given by hkp/m, but in the interacting case it is renormalized by interactions between the (dynamical) particles 
close to ±kp and the particles trapped in the Fermi sea. In the Hamiltonian Q the operators are normal ordered 
with respect to the non-interacting Fermi sea. 

Although the k quantum number is in the low energy approximation restricted to small deviations from ±kp, this 
restriction can be lifted. This is so since the low energy sector of the theory is not affected by the extension of the 
values of k. Without the restriction the model describes two types of fermions, characterized by different values of Xi 
both types with linear dispersion. 

The low energy Hamiltonian can be bosonized, in the standard way, by expressing the Fourier components of the 
charge density operators, for g ^ 0, as boson annihilation and creation operators. 



/ ^ / 27r ^ 

= y l^E^(^«)^x,g> 4 = y |^E^(^^)/'x,-<? (5) 

with 9(q) as the Heaviside step function. The q — components of the charge densities define the conserved fermion 
number and chiral (current) quantum number 

^ = E^x-E^4,fcCx,fc:, >/ = Ex^x-E^4,fcCx.fe (6) 

X fcx X kx 

The bosonized form of the Hamiltonian is [22] 

irh 



^EM 



(7) 
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when modified relative to the fermionic Hamiltonian Q by adding terms that are constant or linear in N. The two 
velocity parameters are 



VN^vp + ^{vm + ¥2(0)) , vj^vp + ^{vm - 1/2(0)) 



(8) 



The low energy sector now corresponds to situations where \q\ is efl^ectively restricted to values much smaller than kp 
and where N and J have only small deviations from their ground state values N = J = 0. 
The bosonized Hamiltonian is diagonalized by a Bogoliubov transformation of the form 



Qq = cosh bq + sluh £^q 

= cosh S^q &J + sinh 6_, 



where ^„ is fixed by the relation 



tanh 2^q = — - 



V2iq) 



Vi{q) + 2nhvF 

In terms of the new bosonic operators the Hamiltonian gets the diagonal form 



H = Y, f^^i blbq + ^{vnN^ + vjJ^) 

g#0 



(9) 



(10) 



(11) 



with the frequency ujq given by 



vf - 



2TTh 



(12) 



If in the low energy approximation the interaction potentials can be approximated by constants, Vi{q) ~ Vi(0) and 
^2(9) ~ 1/2(0), then the bosonic Hamiltonian can be given the field theoretic form 



2 



dx 



(13) 



with u = y/vj^vj and g = \fvjJvN- The fields Q{x) and <I>(a;) are related to the bosonic creation and annihilation 
operators in the following way 



e(.) = e„ + -j-z5:-_= 

^ q^Q V^^Lg\q\ 



2TiL\q 



{bqe'"'' - ble-""") 
sgn{q){bqe''^-~byn 



(14) 



with Oq and are x independent operators. 

In this formulation O and dx^ — d^/dx (or alternatively <& and 9x0) are regarded as conjugate field variables, 
with the basic field commutator given as 



[e(x), 9x$] = [$(x), 9x9] = -Six - y) 



(15) 



where S{x) should be interpreted as the periodic delta-function on the circle of length L. The field dx^ can be 
identified with the fermion number density p{x) of the original description while dxQ is proportional to the current 
density j{x). 

In the representation (13 1 the N and J dependent part of the Hamiltonian (11) does not appear explicitly, but is 
hidden in the zero mode, which is the non-propagating q — part of the theory. This mode is linked to topological 
properties of the fields, which are reflected in the following quasi-periodic condition 



e{x + L) = e{x) + J, <i>{x + L) = ^{x) - N 



(16) 
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where N and J are integers that are restricted by the fermion condition that N ± J are even. In the field theoretic 
description the fermions can therefore be interpreted as topological excitations of the (<i>, 0) fields, and the N and J 



dependent part of the Hamiltonian (111 can thus be interpreted as a topological term. 

In the expansion of the field operators <I>(x) and Q{x) the N and J dependent terms, but also the constant operators 
$0 and Go, represent the zero mode. The presence of <i>o and Oq is needed in order to satisfy the canonical commutation 



relations (15 1, and from this follows the following commutators 



[eo,iv] = -[a>o,j] = — (17) 

TT 

They show that Qq and $o generate operators that change the fermion numbers. The two fundamental fermion 
creation operators, which increase either N-^- or iV_ by one unit, are expressed in terms of the operators Qq and <f>o 
as exp(i7r(6o + $0)) and exp(i7r(6o — *i>o))- 

III. CHIRAL SEPARATION 

The fields Q{x) and ^(x) satisfy a one-dimensional wave equation, and can be separated in a natural way in two 
parts, its right- and left-moving components. These components, which can be defined as 

e±{x) = Q{x)T-Hx) (18) 
9 

satisfy the linear differential equations 

idt±ud,)e±{x) = o (19) 

The (Fourier) expansion of the chiral fields have the following form 



e±(x) - eo± ± ^Q± - 1 y ^z^(^?^''' - ^^'^^ (20) 



with the zero mode operators 



eo±-eoT-$o, Q±^liN±gJ) (21) 
9 2 

These operators satisfy the ^-independent commutation relations 

[QQ^,QA^--S^J (22) 
TT 

with i, j — ±. 

The operators Q± have a natural interpretation as the chiral fermion number operators. For g 1 they take 
non-integer values and for g = 1 they coincide with the operators N±, which take integer values. However, as we shall 
discuss further there are some complications concerning this interpretation of the operators. 

It is then of interest to have a closer look at the topological sector of the theory, which is described by the zero 
mode operators. The quasi-periodic conditions satisfied by Q{x) and <&(a:), can in a natural way be interpreted as 
the condition that the variables Qq and <i>o define a two-dimensional space with the topology of a torus. The precise 
periodicity condition is given by the identification 

(eo,$o) = (eo + n-m,4>o + " + TO) (23) 
with n and m as (independent) integers, and only functions of Qq and <f>o that respect this periodicity are to be 



considered as observables. The periodicity (23 1 is dictated by the spectrum of the two operators N and J, and it can 
be interpreted as a condition on the fermion creation operators exp(i(Oo ± 'i'o)); demanding that they preserve this 
spectrum. 

As an important point to note, the distinction between the two types of independent fermion creation operators 
matches the separation into the two types of chiralities when g — 1. Thus the operators exp(z(Oo ± <i>o)) create 
fermions with well defined chirality. However, for g ^ 1 that is no longer the case, and the corresponding mismatch 
can be seen as a confiict between the dynamical separation of the two chiralities and the periodicity requirement. The 
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former determines the form of the operators Q±{x) associated with the right- and left-going modes, while the latter 
determines the form of the operators exp(j(Oo±<l>o)) that produce changes in the fermion numbers. As a consequence 
of this the operators exp(«0o±), which at the formal level create states with fractional fermion numbers Q±, are not 
acceptable as creation operators since they do not respect the periodicity constraints that the physical Hilbert space 
states should satisfy. Only certain combinations of these operators satisfy the periodicity conditions and therefore 
create states with acceptable (integer) fermion numbers. In Ref. [15] it is argued that operators can nevertheless be 
defined that create chiral excitations with fractional charge and statistics. However, these are well-defined operators 
only if the constraints introduced by the boundary conditions are lifted, and conclusions based on the representation 
of these operators are therefore not fully convincing. 

For the further discussion we find it useful to associate two different quantum numbers x ^^'^ -T with these two 
sides of the chiral separation. In the original fermion model x characterizes the two branches of fermion states, and 
is therefore linked to the fermion numbers N-^ which are fixed by the boundary conditions. In the non-interacting 
case X ^Iso defines the chirality of the state, so that x = -1-1 corresponds to right-moving and X = — 1 corresponds to 
left-moving modes. When interaction is turned on this picture changes. The fermion numbers are still conserved, 
but X is no longer a true chirality quantum number, since a right-moving mode will not be a pure x = +1 mode, but 
also involve a (small) component of x = ~1- This mixing of the two fermion branches is caused by the interaction, 
and is represented by the Bogoliubov transformation which diagonalizes the bosonic Hamiltonian. For the interacting 
system we therefore specify the chiral modes instead by F, so that F = -I- corresponds to the right-going and F = — 
to the left-going modes. Obviously for g — I we have x = T. 

In order to illustrate the meaning of the two quantum numbers, we consider a two-dimensional representation 
of the many-fermion system. As discussed in Refs. ESj EH a two dimensional electron gas in a strong magnetic 
field, when confined to the lowest Landau level, is equivalent to a one-dimensional system, and with a harmonic 
confinement potential in the direction orthogonal to the a;-axis (j/-direction) the Hamiltonian gets the form ([T|) when 
mapped to one dimension. In this case the interaction, in its one-dimensional form, is non-local and therefore gives 
Vi (x) ^ V2 (x) . The same correspondence has recently also been applied in Ref. [T8j where the possibility of detecting 
charge fractionalization on the edge of a quantum Hall system has been analyzed. 

The mapping between the densities in the one and two dimensional representations is given by 

P2{x,y)=Af J dx' J d^e-^[(---')'+i«'-^«lpi(^'+f,^'-f) (24) 

where A/" is a normalization factor, I is the magnetic length of the electron system and pi is the off-diagonal density 
operator of the one-dimensional system. The above expression shows that p2{x,y) is closely related to the Wigner 
function W{x^ k) of the one-dimensional system, 

P2{x,y) ^N' j dx' j dA:e-i^[("-"')'+(''-'='')'lW^(x',fc) (25) 

However P2{x^ y) is, as opposed to W{x, k), a positive (semi-definite) operator for arbitrary (x, y) due to the Gaussian 



factor in (25 1, and has therefore the character of a true density function. The expression also shows that the y- 
direction in the plane is essentially the momentum direction in the phase space of the one-dimensional system, with 
the identification y — kP. The two-dimensional description can thus be considered a phase space representation of 
the one-dimensional system, with a particular definition of the two-dimensional density, and with the Landau-level 
model as a specific, physical realization. 

The filled Fermi sea, in this representation corresponds to a band of integer Landau level filling between two parallel 
edges, located at y = zLkpP = iyp- For the non- interacting system the two-dimensional particle density is explicitly 
given by the x-independent function 



y+ kpf 



erf 



y 



P 



(26) 



with erf(z) denoting the error function. When the interaction is turned on, the ground-state density is essentially 
given by the same function, but with a slight softening of the edges at ±yp [2^ . 

If an additional particle is now injected into the system at the upper (lower) edge that will introduce a density 
modulation of the edge that, for the non-interacting system, will travel to the right (left). For the interacting system 
a density modulation introduced on the upper edge will instead separate into both a right-moving and a (smaller) 
left-moving component. Since the charge on each edge is separately conserved, the separation of the modulation into 
a right- and left-going mode has to be accompanied by a further splitting of each of these into components on both 
edges. The situation is illustrated qualitatively in Fig.l, where an initial wave packet of Gaussian form on the upper 
edge (x = +) dynamically splits into a dominant right-moving packet and a smaller left-moving one. 



FIG. 1: Two-dimensional representation of chiral separation in an interacting system. A charge which is injected at the upper 
edge (with momentum close to kp, blue dashed curve) is separated into a a right-moving and a left-moving component. Each 
of these carries a density component both on the upper and on the lower edge. The ratio between the upper and lower part 
of the right-going charge is determined by the interaction parameter g, and is the same as the ratio between the lower and 
upper part of the left-going charge. The initial condition restricts the total charge on the lower edge to be zero. The figure is 
qualitative and for illustration purpose the interactions have been chosen rather strong with a parameter value about g = 0.5. 



The two-dimensional representation indicates that the injected charge is separated in four parts, which we may 
specify by the two quantum numbers x ^^'^ T. The first one (x = ±) specifies the location on either the upper edge (at 
or the lower edge (at —kp), while the second quantum number F = ± distinguishes between the right-moving 
and left- moving components. For a well localized initial charge distribution, like the one illustrated in Fig.l these 
components are all well defined. There is thus a dynamical separation of the charge in the right- and left-going parts, 
and a frequency separation can further distinguish the parts that are close to +kF and —kp. However, one should note 
that for a general charge distribution the total (integrated) charge has a clear separation only in the two conserved 
charges N-^ corresponding to particles on the upper/lower edge. There is in general no unique definition of the total 
right-moving or left-moving charge. This is an important point that we will now discuss further. 

We first note that in the Fourier decomposition of the charge density the information about the total (global) 
charge sits in the q — Q component, while the local distribution of charge is determined by the q ^ components. 



The periodicity conditions (23 1, which restrict the charges of the upper and lower edges to integer values, only 
affect the q = component. On the other hand the separation in right- and left-going components is only meaningful 
for q ^ 0, since the ^-independent q — component has no motion, neither to the right nor to the left. Therefore, 



the separation of the total charge N into the two chiral components Q±, as done in (21 1, is in this sense somewhat 
arbitrary. In the bosonization of the theory the q = and the q ^ components are treated in different ways, and 
the main point in the following will be to do that also in the discussion of the chiral separation of charge. 

Following Ref flQ] we then introduce a distinction between the total local charge, defined as the g — > component of 
the charge density, and the global charge, defined as the g = component. For the two charges associated with the 
quantum number x the local charges are denoted 

^x = l™,px, (27) 

marked with a bar to indicate that they are not necessary equal to the global charges — Pq=o. Note that the local 
charges, when defined in this way, do not refer to localization within a fixed line element Ax. They should rather be 
viewed as the charge operators restricted to a line element Ax, in the limit Ax — > oo, when this is taken after the 



limit L — > oo. (A more specific way of doing this will be discussed later.) We shall assume that the limit (27) does 
exist for the states of interest. As opposed to the global charges N^, the local charges will take values that are 
not restricted by the boundary conditions. The difference between the two types of charges is due to the presence 
of a compensating background charge that affects only the q — mode. In the case of a compact space, which we 
consider here, this background charge will be evenly distributed over the circle. For an open space this compensating 
charge can instead be viewed as a boundary charge infinitely far away from the region of interest [T^l . 

To find explicit expressions for the local chiral charges, we return to the expressions ([s]) for bosonic annihilation 
and creation operators, extract the charge density operators and re-express them in terms of the transformed Bose 
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operators bq and 6j (defined in (j9|). For q ^ the operators are 



^ {e{xq) cosh + e{^xq) sinh Q \ 



^{e{xq) sinhC, + cosh^,) 6% (28) 

Since 6^ defines the positive frequency part of the operator, with time evolution e"*'^'*, and 5^ defines the negative 
frequency part, with time evolution e'"«*, the time evolution of the operator p^q (for g 7^ 0) is implicitly given by the 
expression (28 1, and the right-and left-moving parts of the operator can therefore be extracted. 



Pxr ,^\H^ [^(rg) Wxq) cosh + d{-xq) sinh Q 6, 



+ e{-Tq){9{xq) sinhe, + 9{-xq) cosh^,) 6^,1 (29) 



This gives the four- fold separation of the charge density specified by the two quantum numbers x and T. 

The local charges, which are defined by taking the limit g ^ 0, we assume to be real valued, which means that the 
limit is independent of taking the limit q ^ 0'^ or g ^ 0^ . This further means that it is the hermitian part of the 
charge density (29 1 that is of interest. We write it as 

I (Pxr , + Plr 5) - [^(rx) cosh + e{-Tx) sinh C,] Br|,| (30) 

with Bg defined by 



B. 



" 2 



We now take the limit g ^ 0, and define 



B± = lim B„ . (32) 

g->0± 



This gives for the xF components of the local charge 

Qxr = [Oi^x) cosh Co + 0{-Tx) sinh^olSr , (33) 

or written out separately 

0++ = cosh^oB+ , Q_+ = sinh ^0-8+ 

Q+_ ^ sinh £_oB- , Q— = cosh ^o^- , (34) 
where the mixing parameters are determined by the Luttinger liquid parameter g as 

cosh^o = ^r7=^ ' sinh^o = ^-r^ • (35) 

Since B± are time independent operators, all the four components of the local charge are separately conserved. The 
two X components of the local charge are furthermore given as 

^+ = E Q+r = ^[(5 + 1)B+ + {9- l)B-] 

N-^Y. ^ ^[(.9 - 1)^+ + (.9 + 1)5-] (36) 
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and the two chiral components as 

X 

Q-^T. Qx- = VgB- ■ (37) 

X 

From these expressions follow that, if we introduce local charges N = , J — X^x' have the relation 

Q±^l{N±gJ), (38) 



which is of the same form as in the previous expression ( |21| used for the chiral charges. An important difference is, 
however, that N and J are not, unlike A'^ and J, restricted to integer values. Therefore Q± can also take more general 
values. 



Note that even if the total local charge in (34 1 has been separated in four parts, these depend on only two charge 
operators, B±. The operator i?+ is associated with the right-going component in such a way that there is, in the 
two-dimensional representation, a fixed ratio tanh^o = (g ~ ^)/{9 + 1) between the parts of this component on the 
upper and lower edges. Similarly _B_ determines the charge of the left-going component with the same ratio tanh^p 
between the parts on the lower and upper edges. This means that for a right-moving component, the (small) charge 
moving on the lower edge can be regarded as a (reduced) image charge of the larger one moving on the upper edge. 
It is natural to view this as a polarization effect caused by the (long range) interaction V2 which acts between the two 
edges. 

In the above discussion we have assumed that the limit (7 ^ is well defined for the charge operators Bq. We are 
aware of the problem to make this assumption precise. Our intention, however, is to apply the charge operators only 
to states that have well defined limits for the expectation values these operators. When we later discuss the charge 
fluctuations we have to be more careful and we will then be more specific about how to take the limit. 



IV. EXAMPLES 



We will now illustrate the general discussion by three examples. The first one is the case where an electron is 
injected at one of the Fermi points, and the charge of the electron splits into two non-integer parts which travel in 
opposite directions [TS]. The second case corresponds to the situation where an initial state is created by injecting 
the electron at one Fermi point of the non-interacting system, and the interaction is then adiabatically turned on. 
In a Fermi liquid the corresponding procedure would turn an electron into a Landau quasi-particle of the interacting 
system. The electron charge is separated into a local charge and a background charge [inj when the interaction 
is turned on, but in this case there is no counter-propagating charge component created. In the third example a 
local charge is created as a local polarization charge by applying an external field to the system. When the external 
potential is suddenly turned off the charge splits into two components which travel in opposite directions. 

We derive in this section the expectation values of the local charge components in the three cases, and follow this up 
in the next section, where the charge fluctuations for the same three examples are examined. Our discussion is based on 
the low-energy approximation to the interacting fermion system using the bosonized form. The limitations introduced 
by this approximation is assumed not to change results in any essential way, even if the charge fluctuations can not 
be regarded as a fully low-energy phenomenon. To check this, in the appendix we compute the charge fluctuations in 
the non-interacting fermion system without making any low energy approximation, both for the ground state and for 
the polarization charge induced by an external delta-function potential. 



A. Example 1: Sudden injection at a Fermi point 



We flrst consider the situation where an electron is injected on one of the edges, which means that the momentum 
of the particle is restricted to an interval which is close to either -{-kp or —kp- This case is the one illustrated in 
FigjT The situation has been discussed in the paper by Pham et. al. |15) and we will compare our analysis of 
the chiral separation of the injected charge with theirs. The (global) fermion numbers N^^ are in this case sharply 
deflned, for example to be A+ — 1 and A_ = 0. When the particle is injected locally we expect no difference between 
the expectation value of the local and global charges, {N^ = N^, since a local injection process cannot introduce a 
background charge that is evenly distributed around the whole circle. 
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With \G) as the ground state of the interacting many-fermion system, the state after injection is 



l*> - *;|G> (39) 

where is a fermion creation operator that injects that particle at either the upper edge (% = +) or the lower one 
{x — Expressed in terms of the fermion field operator ipxix) it has the form 

^1= I dx(j){x)4,l{x) (40) 



where 4>{x) is the wave function of the injected particle. The local charges axe determined by the expectation value 
of the operator Bq , which we write as 

= t^"*^] + [*„&a*^JG)), (41) 

where in the last expression we have used the fact that hq\G) = 0. 

Since hq and &J are directly related to the fermion density operator, they have simple commutators with the fermion 
field operator 



[hq,i^{{x)\ = J ^^[9{xq) cosh iq-e{-xq)smhCq]e-''^-i;l{x). (42) 



/ 27r 

= J^^[0{xq)cosh^q-e{-xq)smh^^lq (43) 



¥^„= / dxe-*«"0(x)V'Ua:). (44) 



This gives 



where we introduced 



For the expectation value of Bq this gives 

{G\^xBq%\G) = ^-[9ixq)cosh^q ^ ei-xq)smh^q]i{G\^x vj/I^JG) + (G|vl/^, vi/t|G)) (45) 
which determines the expectation value of B± in the limit q 0, 

(Br) - lim (Gl^-^ Bq ^-t |G) = d{xT) cosh Co - Oi-xT) sinh^o ■ (46) 

We here assumed limq_^o± "if^q — and used the normalization condition (Gj^*^ ^J-IG) — 1. For x = + this gives 
for the expectation values of the local charges, 

{Q++) = cosh^ xo , (<?-+) = sinh cosh , 

{Q+-) = -sinh^xo, (0 — ) = -sinh^ocosh^o (47) 

and from this follows 

(iV+) = l, (iV_) = 0, (48) 
which are the same as the corresponding values of the global fermion numbers N± . For the chiral charges we find 

= coshCo(coshCo +sinhCo) =^(1+5) 

{Q-) = -sinheo(cosheo + sinheo) = ^(1 -.9). (49) 

If the fermion is instead injected at the lower edge, with x = the results are the same, if the signs of both x and 
r are changed. 

The results (49 1 are the same as with the definition (21 1 for the (global) chiral charges Q±, previously used in Ref. 
1151 This coincidence can be understood from the way the charge is introduced; no constant background charge is 
created, and the limit g — s- is therefore continuous. In the two next examples this will not be the case. 



11 



y 




X 



FIG. 2: Two-dimensional representation of the situation discussed in Example 2. A fermion is injected in the non-interacting 
system at the upper edge (i.e., with positive chirality). The interaction is then adiabatically turned on to give a modified initial 
state indicated by the dashed blue curves. This is a positive chirality state which moves to the right during the time evolution 
of the interacting system. In this case no negative chirality component is created. However, constant background charges on 
the two edges compensate for the fact that the two local charges N+ and N- are not identical to the conserved charges A''+ = 1 
and A'^- = 0. The edge positions before injecting the particles are indicated by the dashed red lines. 



B. Example 2: Adiabatically turning on the interaction 

We next consider the situation where the particle is injected into the non-interacting system and the interaction is 
then adiabatically turned on. The energy eigenstates of the non-interacting system are then turned into eigenstates 
of the interacting system and chiral states of the non-interacting system are turned into "dressed" chiral states of the 
interacting system. The initial state can now be written as 

\^) = U¥^\F) (50) 

where is the ground state of the non-interacting system, ^'^ is the same fermion creation operator as in the 
previous example and U is the unitary transformation from the non-interacting to the interacting system. We note 
that the initial state of the previous example can be written in a similar form, but with a change in the order of the 
operators, |^') = ^l.U\F). 

We evaluate again the expectation value of the operator Bq, 

where we have used the relations 

ag = U^,U, al^U^blU (52) 

The expression has the same form as in the previous example, but now for the non-interacting case, corresponding to 
cosh^q = 1 and sinh^, = 0. This gives 

{F\^^U^B,U¥JF) = ^0{xq){F\^^^l, + ^^,^l\F) (53) 

with the following value in the limit (7 — s- 0* , 

(B±) = ^lim (Fl^-x U^BgU^llF) - 9{±x) (54) 

For X — + local charges now are 



(<3++) = cosh^o {B+) = cosh^o , {Q-+) ^ sinh^o {B+) = sinh^o 

(Q+_) =sinheo(S_) =0, (Q__) =cosheo(S_) =0 (55) 
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which further gives 

(7V+)=cosheo = |^, (7V_)=sinh^o= 1^ (56) 

and 

= cosh^o + sinh^o = ^5 
{Q-) = (57) 

We note that the adiabatic switching on of the interaction has transformed the state with positive chirahty of the 
non-interacting system into a state with positive chirahty of the interacting system. So in this case there is no sphtting 
of the initial charge distribution into a right-moving and a left-moving component. Instead a constant background 
charge density has been created to account for the fact that local charge is not preserved under the adiabatic process. 
This gives the inequalities (-/V_|_) ^ iV+ and (-/V.) ^ N_ which correspond to a discontinuous transition q 0. 

The results found here are consistent with those of Ref. H^I where a first order perturbative result is given for the 
local charge of the dressed fermion, with expansion in the deviation .g — 1 from the free theory. 




X 



FIG. 3: Two-dimensional representation of the situation discussed in Example 3. A charge is initially attracted towards the 
point X = by an external potential (dashed blue curve). The size of the local (polarization) charge is determined by the 
strength of the potential. When the potential is next turned off the charge divides in two equal parts with different chiralities 
which move apart in opposite directions. The local charges are compensated by constant background charges of opposite signs. 
The edge positions before injecting the particles are in the figure indicated by the dashed red lines. 



C. Example 3: Polarization charge 

As a third example we consider an initial state which is created by an external potential W{x) that is slowly turned 
on. The potential will attract a local charge. We assume in this case that no external charge is added, so that 
— N_ = 0. The initial state is the ground state of the Hamiltonian with the additional term, 

AH^ fdx Wix) pix) ^^Y^W.p.g^ ^NWo + ^ E ^« Px-g ' (^8) 
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where Wq are the Fourier components of the potential. The first term gives no contribution since = and the 
second term we write in bosonic form as 



AH ^ \Y.^qP-q 

7 E vl^(cosh + sinh Q{W; \ + Wqb\) 



q^a 



with 



2ttL 



(cosh^5 + sinh^q)Wg 



(59) 



(60) 



The total (bosonized) Hamiltonian, in the subspace with N ^ J = 0, then takes the form 



-){bq 



(61) 



where the last term is a constant that can be ignored. The first term has a form which is identical to the Hamiltonian 
without the external potential, except for a constant shift in the operators bq and This shift can be expressed as 
a unitary transformation 

A ... A* 

(62) 



with 



S = exp 



A* A 



(63) 



This gives H' = SHS^ when the constant term in (61 1 is omitted, and the ground state of the Hamiltonian H' , with 
the external potential included, is therefore 



l^*) = exp 



V(^6 -^6t) 



\G). 



(64) 



If the external potential is suddenly turned off, this state becomes the initial state which is then dynamically separated 
in a right-moving and a left-moving component. The situation is illustrated in Fig|3] 

To find the local charges associated with the right-moving and left-moving components, we analyze the system in 
the same way as for the two previous examples. We find 



mBq\^) ^ iyM(G|5t(6, + &t)5|G) 



1 /^ A, + a; 

2 V 27r kuj„ 



■ (cosh S^q + sinh ) Re Wq 



(65) 



where we have used bq\G) = 0. This gives for the limit q ^ 0, 



(66) 
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with Wq = J dxW{x) and u = y/vjVN as the velocity of the chiral modes, 
separates in chiral components of equal value 



{Q+)^{Q-) 



Wo 
2Trhu 



9 



We note that the polarization charge now 



(67) 



The size of the charges thus depend on the strength of the external potential as well as on the parameters of the 
one-dimensional system. As a check on our calculation, we have verified that the result for the total charge found 
here is consistent with the expression for the full polarization tensor of the Tomonaga-Luttinger model which was 
obtained using perturbative methods |25|. This also suggests that the phenomenon of chiral separation more generally 
could be understood as a dynamical polarization effect. 

Also the local charges associated with the two values of the parameter x have equal values, 



(^±) = {Q±) 

and they are different from the global charges N± — 0. 



Wo 
2'Khu 
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(68) 



CHARGE FLUCTUATIONS IN THE ONE-DIMENSIONAL SYSTEM 



An interesting and important question is to what extent the separation of a local charge in a right-moving and a 
left-moving component should be viewed as a separation of the charge itself or rather as a splitting of the probability 
amplitude, so that the charge moves either to the right or left without being split in two. The answer to this is not 
obvious when the local charge sits on the top of the background charge distribution of the many-particle ground state. 
In such a situation the charge fluctuations of the ground state tend to mask the sharpness of the additional local 
charge. 

It is known from earlier discussions of the phenomenon of charge fractionalization in one-dimensional systems 
[T2t [T3l [26] that in some cases the effect of the background fluctuations can be filtered out by introducing a soft 
sampling function in the definition of the local charge. This happens when there is a gap in the energy spectrum so 
that low frequency fiuctuations are effectively suppressed by the gap while the soft edges of the sampling function 
makes it insensitive to high frequency fluctuations. The sharpness is then defined by the variance of the local charge 
collected by the sampling function over a finite region of the one-dimensional space. 

In the present case the situation is different, since there is no gap in the low-energy spectrum. The background 
fluctuations therefore cannot be suppressed completely by defining a sufficiently soft sampling function. However, 
when sampling over a finite region the variance of the ground-state charge can be reduced to a finite value, and the 
variance of the additional local charge of the excited state can be measured relative to this value. In fact, for a system 
with a finite Fermi momentum kp the ground state fluctuations are flnite (but large) even for a sampling function 
with sharp boundaries, and the effect of these can therefore in principle be subtracted. Our approach will therefore 
be to measure the fluctuations of the local charge relative to the ground-state fluctuations. 

The local charges, which we examine here, are always affected by the background fluctuations, while the global 
charges are not. Since the (sharp) values of the global charges are fixed by the boundary conditions, they are 
insensitive to the dynamical splitting of the charge, and also to polarization effects. Only for modified boundary 
conditions can the sharp charge values therefore be fractional. The important point is that the local charges are not 
restricted in any similar way. 

In the discussion to follow we first introduce a specific sampling function, which collects the local charge over a 
finite interval a. We then use this to evaluate the ground-state fiuctuations of the system and show that these diverge 
both in the infrared and in the ultraviolet. We next examine the charge fiuctuations in the three examples discussed 
previously and show that for sufficiently large a the fluctuations are identical to those of the ground state. As the last 
part of this section we relate the sharpness of the fractional charges in these cases to the fact that the initial state, in 
all the three cases, can be viewed as a coherent state of the bosonic variables in the long wave-length limit a oo. 
For the case of the polarization charge, this is the case even for finite values of a. 

We choose the following form for the sampling function 



/(x;a, b) 



1 



erf 



X + a/2 
b 



— erf 



X — a/2 
b 



(69) 



which essentially equals 1 for —a/2 < x < a/2, and has an exponentially fast fall-off to outside this region, with b 
as the characteristic length scale of the transition. The function f{x; a, b) is well defined on the circle when a << L, 
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and has the Fourier transform 



fg{a, b) = / dxfix; a, 5)e-*«^ = - sm{aq/2)e-'''''' . (70) 



The corresponding charge operator, for a given given chirality T, is 

Qr{a,b) = dxf{x;a,b)pr{x) 



L 

4 ^ oiT,f^^^e-h^^\cosH, + sinh^,) . (71) 



In this expression we have disregarded the contribution from the 5 = term, which does not have a well defined 
separation in the two chiralities. In the limit L — > cx) with a fixed, this term vanishes, and the expression for the 
charge operator is 

Qr(a,6)=^ r e-i^^''^(coshe, + sinhg,) Br, . (72) 

Jo q 

We note that the expression earlier introduced for the (total) local charge Qr is now recovered if, as the next step, 
we take the limit a ^ oo with b fixed, assuming a smooth transition lim^^o± Bq = B±, 

lim Qria, b) = (cosh^o + sinh^o) Br = Or (73) 

a — *oo 

This follows since large a means that the integral only gets contributions from small q, with the following limit for 
the integral 

f°° , sin(ag/2) _i„2j2 tt 
hm / dq — ^ ' ' e 49 " ^ _ (74) 
a^oo Jo q 2 

However, for the variance of the charge operator the limit a ^ 00 is not necessarily a smooth limit. In particular 
will the variance of the ground state charge diverge, and we examine therefore these fiuctuations for finite a. For the 
excited states the charge fiuctuations are finite when the ground-state contribution is subtracted, and with this in 
mind we shall still use the expression Qr = hma^oo Qr{a, b) for the total local charge. 



Ground-state fluctuations 



For the ground state the expectation value of the charge operator (72 1 vanishes since bq\G) =0 and therefore 



{G\Bq\G) = yM^G\bq + bl\G)=0. (75) 



The expression for the variance is, in the limit L — > 00, 

AQr{a,bf = {G\Ql{a,b)\G) 

2 



-^^ Jo q 



« 4 rdr;^e-V(Va)^ (76) 
Jo V 

where in the last step we have assumed a to be sufficiently large so that £,q can be replaced by its g = value, with 
g = (cosh^o + sinh^o)^- The integral can be solved in terms of a generalized hypergeometric function and gives 

AQria, bf ^ 2^2(11, 1}, {3/2, 2}, -a' /2b') (77) 
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FIG. 4: Ground state fluctuations of the charge operators Qr{a,b). The plot shows the variance A.Qr(a,b)^ as a function of 
a/b with the normahzation g — 1 for the interaction parameter. 

This result shows that AQr{a, b)^ is independent of the chirality F = ±, as it should be, due to left-right symmetry. 
The fluctuations only depend on the ratio a/b, as follows from the scale invariance of the gapless theory, and the 
functional form is shown in Fig|4] Note that the interaction parameter g only appears as a normalization factor. For 
large a/b the function has a logarithmic behavior 

AQria,br^^H^). (78) 

The expression shows the expected logarithmic ultraviolet divergence when 6 — > 0, which corresponds to the situation 
where the sampling function approaches a step function. However, there is also logarithmic infrared divergence when 
a —^ oo. In reality the ultraviolet behavior is an artifact of the approximation used in the boson representation. 
When the finite depth of the Fermi sea is taken into account that will introduce a physical cutoff, namely the Fermi 
momentum kp- This cutoff essentially replaces b with l/kp in the logarithm when — )■ 0. The infrared divergence 
for large a is however real, and is due to the presence of massless excitations in the system. 

B. Charge separation and charge fluctuations 

We consider next the three examples previously discussed, where a local charge is dynamically separated in a right- 
going and a left-going component. In all three cases we examine the fluctuations of the total local charge Qr = \/gBr 
of each of the chiral components. Due to the presence of ground-state fluctuations we do not expect the fluctuations 
to vanish, even if the charges of each of the two chiral components in some sense are sharp. We shall here define the 
charges to be sharp if the fluctuations are identical to those of the ground state. 

Sudden injection 

The initial state in this case is j^*) = ^'J,|G) with ^fj, — J dx (t){x)^p^{x), where (f>{x) is the wave function of the injected 
particle and x specifies whether the particle is injected with momentum close to kp or —kp . Since the expectation 
value of the charge operator is already evaluated, we only need to determine the expectation value of the quadratic 
operator, and we examine therefore the expectation value of the operator BqBqi, 



^\\G)\ (79) 
Here we again have used that bq annihilates the ground state. The commutators in this expression can be evaluated 




-(G| 



[fe„vI/t]|G)+(G| 
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by use of ( 43 1 to give 



on 4 



X {(G|^x^l,(,+,olG) + (G|*x.*L'|G) + ^•^•} (80) 



where ^'^ ^ is defined by ( 44 ) . The first term in this expression is identical to the ground-state expectation value of 
the operator. In the limit L ^ oo the expectation value of the squared charge operator then has the form 

(vPlQrK = A rdg ^^"'(°g/^) e-^^^^^(coshg,+sinhe,)^ 

^ Hdq r^^. Bin(ag/2)sin(agV2) ^^_.,.(,.^,., 



q q' 

X (6'(xr) cosh cosh + ^(-xr) sinh sinh^,/) 
X (cosh + sinh ^g) (cosh S,q' + sinh S^qi ) 

X {(G|*x*I(r,w)l^) + (G|*x(r,)*i(r,olG) + ^■^] 

One should note the different behavior of the ground state contribution, on the first line, and that of the remaining 
contribution when the limit a — s- cx) is taken. In the first term the oscillating factor does not give an effective limitation 
in the contribution to the q-integral for large q. The integral is instead limited by the exponential factor, so that 
effectively one has g < 1/h. Furthermore the 1/q factor gives rise to the logarithmic dependence of a/6. The oscillating 
factors in the second term instead introduces the effective limitations g, <?' ^ 1/a. For the integrand the large a limit 
can therefore be interpreted as g' — > 0, and this gives a finite contribution when a — > cx3. The expectation value in 
this limit therefore simplifies to 

(Or) = (O^)o + (0(xr)cosh2^o + ^(-xr)sinh2^o)(cosh^o+sinh^o)' 

= {Ql), + {Qrf (82) 

where the term labeled is the divergent ground state contribution. This expression shows that the charge fluctuations 
of the state \^) = ^'J.jG'), in the limit oo, are identical to those of the ground state, 

AQ^ = (AQj,)2 . (83) 

This result we expect to be true not only when a tend to infinity but when a is much larger than the width of the 
wave function 4>{x) of the injected particle. 

We conclude that when a fermion is injected with a sharp value for Xj it will spontaneously split in two parts of 
different chirality and with non-integer charges ^(1 ±5), where each of the charges is sharp in the meaning that the 
charge fluctuations are identical to those in the ground state. 

Adiabatic initial state 

The initial state in this case is j^*) — t/^'^ji^) with \F) as the non-interacting ground state and U as the unitary 
transformation to the interacting system. As discussed in Sect. 4. 2, the effect of this transformation is to replace the 
matrix elements of Bq for the interacting system with the corresponding ones for the non-interacting system. In the 
present case this means that ( 80 1 is replaced by 



{F\^^u^BqBq,u^{\F) = -Ms^^, + -e{xq)e{xq') 



X 



+ (Pl'^x.^lq'lF) + c.c] (84) 
The corresponding expression for the expectation value of (5p(a, b) gives in the limit a — > 00 

(Qr) = (Q^)o + e(xr)(cosh^o + sinhCo)' 

- {Qt)o + {Qt)' (85) 
from which the same equality follows as in the previous example, 

AQl = (AQr)2 (86) 
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In the present case the initial condition gives rise to a propagating wave with only one chirality, which carries the 
charge It is also here sharp in the sense that the fluctuations are identical to those in the ground state. 

Polarization charge 

This is the case where a local charge is created not by injecting a particle from the outside but by polarizing the 
one-dimensional system with an external potential. The initial state is now j^*) = exp '^„ij-j-bq — -^r^J) \G), with 



the coeflicient Aq defined by (60 1. The state is in fact a coherent state, since the effect of the unitary transformation 
that maps from \G) to j^*) is simply to introduce a c-number addition to the operators bq and as shown by (62 1. 
This means that for any linear combination of b and b^ operators the fluctuations are the same as in the ground state. 
This is in particular so for the operators Bq and for a product of two operators we have the identity 

MBqBq,\M>) - (Vl/Ii3,|vl/)(V|/|B^,|V|/) ^ {G\BqBq,\G) - {G\B q\G) {G\B q, \G) (87) 

For the charge fluctuations this implies 

AQl{a,b)^{AQl{a,b))o 



The conclusion is that also here the fractional charges associated with the two chiral components are sharp in the 
sense that the fluctuations are identical to those in the ground state. However, one should note that the fluctuations 
in this case are identical with those in the ground state, not only for large a, but for any a and b, even if this means 
that the sampling function catches only a fraction of the full charge distribution. 

It is of interest to note that, as opposed to the previous two examples, the chiral charges {Q±) in this case may be 
fractional even in the non- interacting case (i.e., for g = 1). This is so, since the value of the charges depends not only 



on g but also on the strength of the external potential, as shown by ( 68 1 . That chiral separation of sharp fractional 
charges is not an exclusive property of the Luttinger liquid, but occurs also in the free Fermi gas, emphasizes that it 
naturally can be seen as a polarization phenomenon. In the appendix we make use of this by examining the polarization 
effect for non-interacting fermions without applying the low energy approximation of the boson representation of the 
theory. A comparison of the results here and of those in the appendix is of interest in order to check explicitly that 
the high frequency contributions to the fluctuations do not in any essential way change the results of this section. 



C. Coherent states and fluctuations 



The sharpness of the polarization charges discussed above can be ascribed to that the state created by the external 
potential is, in the low energy approximation, a coherent state in the boson variables. It is in fact a coherent state 
in all the bosonic modes labeled by the frequency variable q, and this is the reason that not only the integrated local 
charge but even the charge density is a sharp quantum observable. However, also in the two other cases we have 
examined the initial state can, in an approximate sense, be regarded as a coherent state. This is not the case for any q 
component, but it is true in the limit q 0, which can be seen as a reason for the sharpness of the total local charge, 
while in these cases the charge density is not sharp in the same sense. We shall discuss this point a bit further, first 
for the case of a sudden injection of a fermion at one of the Fermi points. 

We consider then the action of bq on the initial state = ^']^|G), 

bq¥^\G)=[bq,^ijl]\G) (89) 



where bq annihilates the ground state \G). The commutator has previously been given in Eq.(42l, and in the limit 
q ^ this gives 



with 



bq^UG) ^ P±^^l\G) , g^0± (90) 



/ 27r 

'-^[0(x)cosheo-e(-x)sinhCo] (91) 



Eq. (90 1 shows that the initial state |^) is an eigenstate of bq, and therefore (approximately) a coherent state, for 
sufficiently small q. This implies that observables that are linear combinations of bq and 6j operators for sufficiently 
small q, will have quantum fluctuations in this state that are equal to those in the ground state. This is in particular 
true for the total local charge operator Qr- 
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In the case where the interaction is adiabatically turned on the initial state is l^f) = where \F) is the 

ground state of the non-interacting system and U is the unitary operator that transforms from the Hamiltonian of 
the non-interacting system into the Hamiltonian of the interacting system. In this case we have 

bgU¥jF) = Uag¥jF)^a±^U¥^\G), g ^ 0± (92) 

where the eigenvalue a±-^ is the non-interacting (g — I) version of the coefficient (3±^, 

(93) 

Since the initial state j^*) can be regarded as an eigenstate of the bosonic annihilation operators bq for sufficiently 
small q, the situation is the same as for sudden injection of the fermion. In both cases the form of the initial state as 
a coherent state for small q can be seen as an explanation for the fact that the fluctuations of the total local charges 
are identical to those of the ground state. 

In both cases discussed above, with a sudden or adiabatic injection of an additional particle, the wave function of 
the particle has been assumed to have a well defined value of the quantum number x- This means that the particle 
is injected with momentum close to one of the Fermi points kp or —kp. As a contrast to these situations let us now 
consider the case where a fermion is instead injected as a superposition of components with x = + and x = — . This 
corresponds to replace, in the above expressions, the fermion creation operator with a composite operator of the 
form = X^y'^x^y- "^^"^ corresponding states ^'^|G) and U'i'^F) then are no longer coherent states in the sense 



expressed by ( pjT| and (93). This follows since the coefficients /?y and have different values for x — + x = ^• 
The variance of the chiral charge in both cases can be written as 

AQl = (AQr)g + Ed^xP " Wxl') (Or)' - 2\c+\'K\^ (Qr)+ (Or)_ (94) 

X 

with (Qr)^ = (G|vI'y(3r*y|G) in the case of sudden injection and (Qr)^ = {Fl"^ Q^^^xl^) '^^ ^^'^ adiabatic case. 
Unless one of the coefficients Cy has absolute value 1 and the other is zero, there will now be non- vanishing deviations 
in the variance from the ground state value, and the mean value of the observable Qr in the state |*) = does 
therefore not correspond to a quantum mechanically sharp value. The chiral separation of the initial charge in this 
case should therefore not be seen as a splitting of the charge itself but rather of the probability for the charge to move 
to the right or to the left. 

In fact, it is clear from general reasoning that any state which has quantum mechanically sharp values for the 
charges Qr (as well as for N^) has to be characterized by well-defined, quantized values for the global charges N^. 
The reason for this is the following. The local charges Qr commute with the global charges A^^. This implies that the 
matrix elements of Qr as well as of Qp vanish between states that have different quantized values for either or 
N^. Consequently, for a state j^*) which is a superposition of such states, the expectation values of Qr and Qr will 
only dep end on the absolute square of the expansion coefficients of j^*) in eigenstates of with different eigenvalues. 



Eq. (|94|) shows this for the particular case where a single particle is added to the system. Also the variance AQp will 
therefore only depend on the absolute square of the expansion coefficients, which means that it has the same form as 
for a statistical mixture of the same eigenstates. This implies that the charges Qr, for such a superposition are not 
quantum mechanically sharp (unless the mean value of the operators are precisely the same for the components with 
different values of N^). 

This makes it possible to illustrate in a simple way the difference between the two situations where, on one hand 
chiral separation gives rise to a separation of the probability for the charge to move in one direction or the other, 
and on the other hand the situation where the charge is split in two sharply defined smaller charges with opposite 
chiralities. We may consider the case of non-interacting fermions, where the first case corresponds to making a sudden 
injection of a particle in a superposition of the two chiralities. The charge is then not split in two, but the particle 
moves with a non- vanishing probability for each of the two chiralities. In the other case a polarization charge is created 
by an external potential. When the charge is released by turning off the potential it divides in two chiral components, 
and both of these are sharp in the meaning that the charge fluctuations are identical to those of the ground state. 
The expectation values (Qr) may be equal in theses two cases, but the fluctuations AQr will be different. 



VI. SUMMARY AND OUTLOOK 



In this paper we used the concept of local charge to discuss the character of the chiral separation of charges that 
takes place in one-dimensional Luttinger liquids. Depending on the initial conditions the expectation values for the 
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chiral components of the charges may take different, non-integer values, and we have examined to what extent these 
charges can be regarded as quantum mechanically sharp. Since the system has a gapless spectrum the question of 
quantum sharpness of an observable is less clear than for systems with a gapped spectrum. We have here defined 
sharp to mean that the quantum fluctuations of an observable in a given state are identical to the fluctuations of the 
observable in the many-particle ground state. 

We have examined the mean values and fluctuations of the chiral charge components with three different initial 
conditions for the quantum state. The first one corresponds to a sudden injection of a fermion at one of the Fermi 
points, the second one to injection of a particle in the non-interacting system followed by an adiabatic switching on of 
the interaction, and the third one corresponds to creation of a local polarization charge by an external potential. In all 
three cases the charges of the two chiral components of the quantum state have been found to have non-integral values 
and to be sharp in the sense that the charge fluctuations are identical to those of the ground state. In particular we 
have found that sharp chiral charges of arbitrary values can be created as polarization charges by applying an external 
potential. This is true also for the non-interacting Fermi gas which shows that chiral separation is a polarization 
phenomenon which is not exclusive to Luttinger liquids. 

The quantum state created by the external potential has been shown to be a coherent state in the bosonic variables. 
Also in the two other cases the quantum state is a coherent state in an approximate sense, restricted to the bosonic 
variables of low frequency. This property, we believe, provides an explanation for the charge fluctuations to be 
identical to the ground state fluctuations. However, since the representation of the states as coherent states is linked 
to the low-frequency bosonized formulation of the theory, this raises the question if the same conclusion concerning 
the sharpness of the non-integer charges can be drawn for the fermion system with a finite Fermi momentum, where 
the high frequency contributions are different. The comparisons we have made in the appendex for the case of non- 
interacting fermions indicate that the high frequency effects will in practice not change the conclusion concerning the 
sharpness of the non-integer charges. 

A related question concerns the importance of the one-dimensionality for the results. The non-integer charges 
studied here can be viewed as arising as due to polarization effects in the many-fermion system. Clearly there 
are related effects also in higher dimensions, and the question is whether non-integer charges which are quantum 
mechanically sharp, in the meaning used here, can be present also there. This problem has been addressed for 
systems with a gap[23], but we know of no analysis of the gapless case. The close link to the coherent states of 
the boson representation may indeed suggest that one-dimensionality is important for the effect, but apart from this 
observation we will have to leave this as an interesting open question. 

Although the examples we have analyzed show that quantum mechanically sharp charges of any value can in 
principle be created, this does not mean that sharp quantum values is a typical feature for the local charges of 
arbitrary quantum states. In the examples we have discussed the quantum states are characterized by quantum 
well-defined, sharp values for the global fermion numbers N^. In particular, in the two first examples, a single fermion 
is added with well defined value of the quantum number x- If this particle is instead injected as a superposition of 
components with both x = + and x = ~ the local charges of the right-and left-going modes will no longer be sharp, 
but will have charge fluctuations of the same form as for a statistical mixture of the two components with different 
values of x- 

In this paper we have not discussed how to actually observe the fractional charges in question, but this most 
important question has been addressed in some recent papers. In the experiment described in Refs. ll6l and fT7l electrons 
are injected by tunneling into a quantum wire at one of the Fermi points, and the chiral charge fractionalization is 
determined from a left-right asymmetry in the resulting current. The experimental result is consistent with the 
prediction of Ref. [1^] and thus also with our result ( 49 1 for the value of the local chiral charges. (We are aware 
that a different conclusion concerning this experiment is reached in Ref. [27l) It is not clear to us that this kind of 
measurement can be used to infer the presence of sharp local charges, but we note that the conclusion in Ref. [T7| 
about sharpness of the charge is the same as reached here, there based on a calculation of the zero frequency noise in 
the current. 

Noise measurements is also the subject of a recent theoretical paper by Berg et al. which considers the detection 
of fractional charges in a physical system very similar to the two dimensional model discussed in this paper [18]. 
The difference is that the width of the Hall bar is assumed to change continuously with a constriction in the middle. 
Away from the constriction the system is essentially free, with g = 1, while g < in the middle region. A proposed 
experiment is to inject electrons on the upper edge via a contact, and measure the absolute value and the noise of 
the reflected current on the lower edge. The result of the calculation is consistent with thinking of the process as a 
reflection at a sharp boundary between a lead with g = 1 and a wire with g < 0. Although the situation is similar to 
the one discussed in some detail in this paper, where a unit charge inserted at the upper edge spontaneously separates 
in a right- moving and a left-moving part, the presence of a boundary makes it different. In particular the value of the 
left-moving charge on the lower edge is not the same in the two cases. It would be of interest to extend our analysis 
to this case, with reflection of a fractional charge at a sharp boundary, and in particular to evaluate the fluctuations 
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in the reflected charge. However, as far as the suggested experiment is concerned it is not obvious to us that reflection 
towards a sharp boundary is the best way to model situation. For a sufiiciently smooth constriction, we think that 
our example 2, i.e., a adiabatic switching of the interaction, might in fact be more appropriate. 

The results found in this paper, that states with arbitrary, sharp values for the local charges of the two chiral 
components can (in principle) be created, seem to cast some doubt on ideas that have been proposed, that quasi- 
particles with specific (fractional) values for the fermion numbers are the natural charge carriers of the Luttinger liquid 
[15]. However, we believe that the idea that fractionally charged quasi-particles have a natural place in the description 
of Luttinger liquids may not be ruled out. In particular, when we consider a system of non-interacting fermions, there 
are low-energy excitations which carry charges with non-integer, sharp values. Even so, the fundamental fermions, 
which carry unit charges, play a special role in the low-energy description, since they define a system of free particles. 

The idea that this picture can be generalized to the interacting case, where the low energy theory can be described 
in terms of non-interacting quasi-particles characterized by fractional charge and fractional statistics is an attractive 
one. We refer the reader to the interesting ideas of Isakov j2S] and Wu ^5] which relate the Luttinger Hquid to 
systems of particles with fractional exclusion statistics. It is also interesting to note that Heinonen and Kohn, in the 
paper we have referred to |19j . assume that the fractional charge, calculated in the situation where the interaction is 
adiabatically turned on, is identical to the (local) charge carried by the Landau quasi-particles, and in that context 
they do not make a distinction between quasi-particles in one and higher dimensions. The questions concerning the 
meaning of such quasi-particles with fractional charge in Luttinger liquids remain as interesting questions for future 
investigation, and this goes beyond the discussion of how to characterize these charges as being quantum mechanically 
sharp or not, which has been the main aim of this paper. 



In this appendix, we examine the charge fluctuations of the non-interacting fermion system, first for the filled Fermi 
sea, i.e., with the system in the ground state, and then in the presence of a delta function potential. In the latter 
case we also calculate the mean value of the induced charge. 



Rather than using the boson representation of the theory we apply the full many-particle description. The aim is 
to see explicitly to what extent high frequency fiuctuations lead to different results in the two approaches. We note 
that, compared to the bosonized description, there is an additional frequency parameter present in the many-particle 
approach, namely the Fermi momentum kp. It provides a high frequency cutoff, and it is of interest to see precisely 
how it appears in the expressions. The high frequency cutoff b, introduced earlier as a smoothness parameter of the 
sampling function f(x; a, b), can now be set to zero giving a sampling function with sharp steps. 



APPENDIX A: FLUCTUATIONS IN THE FREE FERMION SYSTEM 



1. 



Ground state fluctuations 




(Al) 



The corresponding charge operator is 




k,k' 



(A2) 



where the Fourier transform of the sampling function is given by 




(A3) 



The ground state expectation value of the operator is 




(A4) 



the quadratic charge operator is 




(A5) 
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and the ground state expectation value of the operator is 

kj^ kf kp 

(Qi) - E E m-r-i. E E fih-h;ar 

ki— — kf' — — kp ki— — kp ^2— — 

^ kp oo 
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-ik{x — x') 



(A6) 



For the charge fluctuations this gives 



L .l-a/2 J-a/2 sin^(7r3^) 



a/2 



(A7) 



By exploiting the symmetries of the integrand and changing to new variables X — -^{x + x') and ^ = x — x' , the 
expression can be re-written as 



sin^(7riVo|;) 



(A8) 



We take the limit a/L ^ with Nq/L = kp/Tr kept constant, and further introduce a new integration variable by 
= arj and a new parameter P = TraNo/L = akp. This gives 



13- 



dr]{l - rf) 



sin p-q 
if' 



The integral can be expressed in terms of special functions as 

1 



AQ^ = — [l + 7 + 7r/3~cos2/3-Ci(2/3)-2/3Si(2/3) + ln2/3] 

— — [1 + 7 + TTakp — cos 2akp — Ci{2akp) — 2akp Si{2akp) + In 2akp] 



(A9) 



(AlO) 



where 7 denotes Euler's constant, Ci the cosine integral function and Si the sine integral function. 




FIG. 5: The ground state fluctuations AQ^ evaluated in the many-ferion formalism for non-interacting particles, shown here as 
a function of akp (solid blue curve). The corresponding result derived by use of the boson representation is also shown (dashed 
red curve). The ultraviolet cutoff parameters in the two cases are identified by the relation 1/6 = 2.5A;f. 
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The corresponding expression in the boson representation is given by (77), with the asymptotic form, for large a/6, 
given by (78 1. In the present case, with 5=1 and with summation over the two chirahties, the asymptotic expression 
is 



AQ(a,6)^ 



For the expression (AlOl the corresponding asymptotic form for large akp is 



1 



lii(2akp) 



(All) 



(A12) 



This shows that the expressions for the charge fluctuations, evaluated in these two ways, agree up to a constant if 
we identify the cutoff parameter b in bosonized case with the inverse of the Fermi momentum in the other case. The 
constant can in fact be absorbed in a relative scale factor for the two cutoff parameters. This is illustrated in the 
figure where the charge fluctuations are plotted as functions of /3 = akp with the identiflcation 1/6 = 2.5 kp. In spite 
of the fact that the high frequency contributions to the fluctuations are not identical in the two cases, the resulting 
curves are very close for all values of (3. 



2. Polarization charge and fluctuations for a delta function potential 

We now turn to the polarization charge induced by an external potential which we for simplicity shall take to be a 
delta function. We examine both the expectation value and fluctuations of the local charge and compare the results 
to that of the low-energy bosonized description. 

The Schrodinger equation reads 

- ^T^^(^) - ^'^(^ - Lm^bix) = E^j{x). (A13) 

where rj measures the strength of the delta function potential located at the point x — L/2, and we take 77 > 
corresponding to an attractive potential. The delta-function potential gives rise to a discontinuity in the derivative of 
the wave function at x = L/2, and it is convenient to restrict the coordinate x of the periodic variable to the interval 
—L/2 < X < L/2, so that this discontinuity simply amounts to imposing the boundary conditions 

f(L/2)-f(-L/2).2,|l*(L/2). (A14) 

and ij{-L/2) =iP{L/2). 

Without the potential (ry = 0) the energy spectrum has a set of doubly degenerate eigenstates, which may be chosen 
as symmetric (even) and antisymmetric (odd) functions, 

i'nix) ^ cos knX , (x) = -y/^ siu fc„a; , n=l,2,... (A15) 

with wave numbers fc„ = 27m/L. In addition there is a single non-degenerate (zero energy) state, which has k — 
and is therefore constant over the circle 

Mx) = ■ (AI6) 

V L 

When the interaction is turned on the zero energy state is changed to a negative energy state. The potential thus 
attracts the charge of the particle in this state and forms a bound state of the form 

ipo{x) = Aocosh{Kx) , (A-17) 

with K given as solution of the transcendental equation 

, fnL 



= rycoth ( — ) , (A18) 
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and the normalization factor given by Aq — y KL+s?nh(KL) • ^^^o the other even eigenvalue functions are modified by 
the potential. The form is the same as without the potential, 



but the values of the momentum variables are shifted and now are solutions of the equation 



m 



-7] cot 



(A19) 



(A20) 



as follows from the boundary condition (A14l. There is a sequence of solutions to this equation, restricted by 

2?™/^ < fc„ < 27r(n+ , n = 0,l,2,... (A21) 
and the modified normalization constants are 
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2 mrjL 
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(A22) 



The odd functions, on the other hand, are not modified by the potential since they all vanish at the point x ~ L/2. 
The degeneracy of the excited states are thus lifted by the potential. 

We consider now the many-particle system where an odd number of particles N = 2M + 1 in the ground state 
occupy the N lowest energy eigenstates, 



M 



iG)=n«^<^)4io) 



(A23) 



with cL as particle creation operators. The field operator is expanded as 



(A24) 



n=l 



T! = l 



and the particle number density is p{x) — {x)ip{x). 

The ground state expectation value of the particle number density is 



k<kF 
M 



M 



(A25) 



71—1 n—l 

and the density-density correlation function 

C{x,y) = {{p{x)~{p{x))){p{y)-{p{y)))) 

= (5(x-y) E i'*ki.x)i^k{x) - E ^ki.x)MxWii.y)i'k{y) ■ 



(A26) 



k<kf 



k.Kkf 



(Note that the background charge of the unperturbed ground state has not been subtracted.) 

We use these expressions to calculate the expectation value and the variance of the local charge operator restricted 
to the region Da of width a centered at x = L/2. To find these quantities we need to evaluate numerically the 
following coefficients 



L/2 

(L-a)/2 
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+ (1 - 5n,n) 



dxr:{^)rm{^) 
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rL/2 

J{L-a)/2 



A* A 



sin((fc„ + k„,)L/2)) - sin((fc„ + fc„i)(L - a)/2) 



+ (1 - ^nm) ( sin((fc„ - fc„0L/2)) - sin((fc„ - ~k^){L - a)/2) 

+ (5„™A;A„^ (A28) 
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A;„i(cosh(KL/2) sm{kmL/2) — cosh(K(L — a) /a) sin(fc„i(i — a) /2)) 

+ K(sinh(KL/2)cos(fc„L/2)-sinh(K(L-a)/2)cos(fc„(L-a)/2)) [■ (A29) 



F = 2 dx'ipQ{x)tpo{x) — — {sinh(KL) — sinh(K(i — a)) + Kfl} . (A30) 

J(L-a)/2 2k 



The expectation value of the fermion number Na is given by 

M M 



(Na) = J2 ^nn + J] + F , (A31) 



n— 1 n— 1 

and the variance is 

M M M M M 

AAT^ — C„„ + D„„ + F — i^^ — CnmCmn + DnrnDrmi — 2 |G„ 

n— 1 n— 1 n,m— 1 n,m— 1 rn — 1 



(A32) 



In order to make a comparison with the corresponding expressions found in the bosonized formulation we include 
two modifications. The first one is to make a subtraction of the constant background charge of the non- interacting 
system and the other is to compensate for the finite value of a/L. The modified expression for the expectation value 
of the local charge is 

1 — a/L 

which is such that in the limit a ^ L, (Ql) = {Nl) = ^• 

In the numerical evaluation lengths are measured in units of L and the strength of the potential in is measured 
by the dimensionless parameter 77 = (mL/}T?)r]. The wave numbers fc„ have been determined by solving numerically 



Eq.(A20l and the value of k has similarly been found by numerically solving (AI8I. These values have been used 
when evaluating the coefficients Dnm, Gm and F. 

The numerically evaluated ground state expectation value of the charge density is shown in Fig. 6 for particle 
number TV = 55 and for two values of the potential strength 77. The polarization charge induced by the delta function 
potential is strongly localized around the point x = L/2, and more so for larger than for the smaller value of rj. The 
high frequency Friedel oscillations in the density occur at the wave vector 2kp = 2'kN /L and their amplitude decrease 
with increasing N . kp also determines the width of the central peak, and this is consistent with the Fermi momentum 
effectively defining an ultraviolet cutoff. 

To compare the results with those found in the boson representation we focus on the value of the total local charge 
Q ~ X)r'3r- In the non-interacting case, with g = 1 and Wq = —rj for the strength of the potential, the value 
previously found for the charge is 
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FIG. 6: Charge density in the background of a delta function potential. The charge distribution is plotted for = 55 particles 
and for two values fj — 1 and ^ = 10 of the potential strength parameter. The density shows in both cases a peak at the 
location of the potential, with largest amplitude for the largest value of fj. The oscillations are a high frequency effect, with 
wave length determined by the Fermi momentum kp- 



We have labeled the charge with "bos" to specify that it is evaluated in the low energy, bosonized approximation. 
The corresponding charge evaluated by using the many-fermion formalism is labeled {Q} f^^m' ^^"^re this charge is 
identified with (Qa) for a value of a that is sufficiently large to capture fully the contribution from the central part of 
the charge distribution. 

In Fig. 7 the charge is shown as function of (Q)^^^- The curve is obtained by varying the potential strength 

77 with fixed particle number, in the figure corresponding to = 201. The curve shows that (0)^g,.„ ^J^d (0)^^^, 
agree well for small values (0)^g^„ ^ 0-1; corresponding to a weak potential strength 77 (measured relative to the 
Fermi velocity). When the strength of the potential increases, however there is an increasing discrepancy. This has 
to be viewed as effects of high frequency contributions, which are treated differently by the two methods. To some 
extent this can be seen as a consequence of using a delta function potential, since this has Fourier components of 
arbitrary high frequency. 

The increase in importance of these high frequency components with the strength of the potential can be understood 
in the boson representation in the following way. The quantum state in the background of the potential is a coherent 
state which can be expanded in powers of the operator that act on the ground state \G). Since is 

proportional to the potential strength r/ this shows that higher powers of and therefore higher energy contributions, 
will be more important for large than for small values of 77. 

The charge fluctuations AQ^ have been evaluated for different values of a/L, 77 and N. The interesting point has 
been to check if, also when evaluated by use of the many-fermion wave function, the fluctuations are identical to those 
of the ground state, so that this result is not just an artifact of the low energy approximation. In Fig. 8 the results 
are plotted as functions of a/L for = 55 and for three different values of 77, one of them corresponding to the value 
77 = of the free system. As shown by the plot there are deviations from this curve for non-vanishing values of fj, 
but these have the character of oscillations about the fluctuations curve of the free system. The oscillations seem to 
be rapidly damped when 77 decreases, and in the plot that is demonstrated by the difference in the amplitude of the 
oscillations for the curves with 77 = 10 and 77 = 100. 
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FIG. 7: Expectation of the total local charge, (Q)^.^^^, determined by the many-fermion calculation, shown as a function of the 

charge {Q)^^^^ determined by calculation in the boson representation. The result is shown by the solid (red) curve for A'' — 201 
particles. The dotted curve shows the line of equality for the two charges. For small charge values (weak potential) the two 
methods to calculate the charge agree well, while for larger charge values (strong potential) they give diverging results. The 
difference can be ascribed to the effect of having a finite value for kp in the fermion calculation. 

It is of interest to compare the results displayed in Fig. 8 with those in Fig. 7, where the results for the evaluations 
of the expectation value of the charge are shown. The discrepancy between the results for (Q)yg^„j E^nd {0})^^^ shown 
in Fig. 7 can be ascribed to the effect of the finite value of the Fermi momentum kp for the fermion result, and the 
curves indicate that this effect decreases rapidly for charge values below (^Q) ~ 0.1. For the particle number TV = 55, 
as used in Fig. 8, the parameter value 77 = 10 corresponds to (Q) « 0.02 and fj = 100 corresponds to (0) « 0.2. The 
curves shown in this figure demonstrate that also the oscillations of the charge fluctuations are rapidly damped below 
(Q) ~ 0.1. It seems natural to assume that the appearance of these oscillations are also due to the finite value of 
kp, and a closer inspection of the periodicities confirms this assumption. The oscillations in the fiuctuations are then 
closely related to the oscillations in the charge density displayed in Fig. 6, and they therefore appear as a consequence 
of the use a sampling function with sharp edges, since a smooth edge would suppress high frequency contributions 
with k K kp. Apart from these oscillations we find no significant difference between the fluctuations of the local 
charge and the fluctuations of the ground state. 

Let us finally stress the point that the problems met here in getting a simpler picture of the charge fluctuations 
are linked to our use of a delta-function potential and a sampling function with sharp edges. Both these functions 
could be smoothened, but that would involve a more demanding calculation than the one we have aimed at in this 
Appendix. 
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N=55 




FIG. 8: Charge fluctuations in tlie background potential. The plot shows the variance AQ^ as a function of a/L for three 
different values of the potential strength r). The particle number is in all three cases A'^ = 55. The deviations from the ground 
state fluctuations are very small for f] = 10, but can be seen in the window with magnification of a section of the curves. 
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